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BUBBLE  DYNAMICS  IN 
A  SUPERHEATED  LIQUID 

by 

W.  T.  Sha  and  V.  L.  Shah 


ABSTRACT 

This  report  presents  an  extensive  literature  survey  on 
bubble  dynamics.  Growth  of  a  single  spherical  bubble  moving  in 
a  uniformly  superheated  liquid  is  considered.  Equations  of 
motion  and  energy  are  presented  in  the  forms  that  take  into 
consideration  the  interaction  between  the  motion  and  the  growth. 
The  fourth-order  Runge-Kutta  method  is  used  to  obtain  a  si¬ 
multaneous  solution  of  equations  of  motion  and  growth  rate, 
and  the  solutionis  compared  with  available  experimental  results. 
Results  for  liquid  sodium  are  presented  for  a  range  of  pres¬ 
sures  and  Jakob  numbers. 


I.  INTRODUCTION 


A.  Bubble  Dynamics 


Broadly  defined,  bubble  dynamics  deals  with  all  phenomena  associated 
with  formation,  growth,  motion,  and  collapse  of  single  or  multiple  bubbles. 
Since  bubbles  are  encountered  in  a  wide  range  of  industrial  equipment,  bubble 
dynamics  is  the  subject  of  study  by  mechanical  engineers,  chemical  engineers, 
geologists,  metallurgists,  and  researchers  whose  main  interest  is  in  boiling 
and  cavitation.  With  the  advent  of  nuclear  reactors,  more  and  more  nuclear 
engineers  are  taking  interest  in  solving  problems  related  to  bubble  dynamics. 

Although  much  work  has  been  done  in  this  field,  the  subject  is  still 
not  completely  understood.  As  a  consequence  of  the  extreme  complexity  of 
the  interaction  between  a  bubble  and  its  surroundings,  most  of  the  design  work 
is  carried  out  on  the  basis  of  experimental  trial-and-error  approaches.  How¬ 
ever,  theoretical  modeling  can  play  a  significant  role  in  providing  physical 
insight  and  rendering  equipment  more  efficient  and  reliable. 

B.  Present  Work 


The  present  study  was  initiated  with  a  desire  of  understanding  the 
dynamics  of  sodium  vapor  bubbles  in  a  fast  breeder  reactor.  After  a  brief 
review  of  past  theoretical  work  relating  to  a  stationary  and  moving  bubble,  a 
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solution  is  presented  for  a  moving  bubble  that  is  simultaneously  growing  in 
a  uniformly  superheated  liquid.  The  fourth-order  Runge-Kutta  method  is 
used  to  obtain  a  simultaneous  solution  of  equations  of  motion  and  growth  rate, 
and  the  solution  is  compared  with  the  available  experimental  results.  Results 
for  liquid  sodium  are  presented  for  a  range  of  pressures  and  Jakob  numbers. 


II.  GROWING  AND  COLLAPSING  PROCESS 
A.  Vapor  Bubbles 


A  bubble  nucleates  in  a  superheated  liquid  due  to  vaporization  of  liquid 
molecules  into  a  space  that  is  unoccupied  by  liquid;  e.g.,  a  cavity  on  the  heat¬ 
ing  wall  or  on  a  suspended  foreign  material.  Once  a  bubble  is  formed,  three 
factors  control  its  further  growth.  These  are  mechanical  forces,  he  at- transfer 
rates,  and  diffusion  rates  at  the  interface.  Initially,  the  growth  is  controlled 
by  inertia  and  surface-tension  forces,  but  as  the  bubble  size  increases,  the 
other  factors  become  more  and  more  important. 

There  are  two  kinds  of  vapor  bubbles:  cavitation  and  boiling.  A  cavi¬ 
tation  bubble  grows  in  a  liquid  that  is  at  a  low  temperature  and  a  negative 
pressure  (i.e.,  pressure  lower  than  a  saturation  pressure).  The  dominating 
factor  controlling  its  growth  is  the  inertia  force.  In  this  case,  a  solution  of 
the  equation  of  motion  along  with  the  continuity  equation  may  be  sufficient  to 
determine  the  growth  or  collapse  rate  of  a  cavitation  bubble. 

A  boiling  bubble,  on  the  other  hand,  grows  in  a  liquid  that  is  at  a  high 
(superheated)  temperature.  In  this  case,  the  heat- transfer  rates  at  the  inter¬ 
face  control  the  growth  rate,  and  the  energy  equation  must  be  included  to 
determine  the  bubble  dynamics  of  a  boiling  bubble.  Moreover,  if  a  bubble  con¬ 
tains  gases  in  addition  to  pure  vapor,  the  determination  must  include  the  dif¬ 
fusion  equation. 

The  study  in  this  report  is  restricted  to  the  growth  of  a  pure  vapor 
boiling  bubble. 

B.  Equilibrium  Conditions 


Consider  a  spherical  vapor  bubble  of  radius  R,  pressure  pv,  and  tem¬ 
perature  Tv,  which  is  in  equilibrium  with  the  surrounding  liquid  of  tempera¬ 
ture  T  and  pressure  p.  (See  Fig.  1.) 

For  thermodynamic  equilibrium  to  exist  at  the  curved  interface,  the 
following  conditions  have  to  be  satisfied. 

1.  Both  the  vapor  and  liquid  have  to  be  in  superheated  states  with  re- 
pect  to  saturation  temperatures  corresponding  to  a  plane  interface;  i.e., 


^  >  ^  satW 


(1) 
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and 


Tv  >  Tsa£(pv).  v2) 

2.  The  vapor  pressure  pv  has  to  be  greater  than  the  liquid  pressure  p 
to  balance  the  additional  surface-tension  force.  The  excess  pressure  in  a 
bubble  of  radius  R  is  given  by  the  Gibbs  equation, 

2a  /  v 

pv-p  =  ir- 

where  a  is  the  surface  tension. 


LIQUID 
Pi  T 


Fig.  1 

Spherical  Bubble  Growing 
in  a  Superheated  Liquid 


C.  Equilibrium  States 


Before  describing  the  thermodynamic 
states  of  a  growing  bubble  on  a  saturation 
curve,  we  shall  consider  the  thermodynamic 
equilibrium  states  of  a  single -component 
fluid.  Figure  Z  shows  a  saturation  curve  on 
a  pressure-volume  diagram  for  a  pure  sub¬ 
stance.  The  curve  ABCDEF  is  an  isotherm 
of  temperature  T.  The  isotherm  can  be  ap¬ 
proximated  by  the  van  der  Waal  equation, 


(Mv)2 


(Mv  -  b) 


RT. 


(4) 


The  curves  AB  and  EF  represent  the 
thermodynamic  states  of  a  single-phase  liquid 


Saturation  Curve  for 
a  Pure  Substance 


and  a  single -phase  vapor,  respectively.  When  both  liquid  and  vapor  are  in 
equilibrium  with  a  plane  interface,  positions  B  and  E  represent  the  states  of 
liquid  and  vapor,  respectively.  The  curves  BC  and  DE  represent  the  unstable 
states  of  superheated  liquid  and  subcooled  vapor,  respectively. 

For  a  spherical  bubble,  the  interface  is  curved.  For  a  curved  inter¬ 
face,  as  described  before,  both  the  liquid  and  vapor  are  in  superheated  states. 
The  state  of  a  surrounding  liquid  at  temperature  T  therefore  lies  along  the 
curve  BC;  the  state  of  a  bubble  vapor  at  temperature  T  lies  on  curve  EF. 


D.  Thermodynamic  States  on  a  Saturation  Curve 

Once  the  equilibrium  conditions  for  a  spherical  bubble  have  been  de¬ 
scribed,  it  is  easy  to  locate  the  states  of  liquid  and  vapor  with  a  curved  inter¬ 
face.  In  Fig.  3,  the  position  L  represents  the  state  of  liquid,  and  the  position  V 

represents  the  state  of  vapor. 


Fig.  3.  Equilibrium  States  of  Liquid  and  Vapor  with  a  Curved  Interface 
It  can  be  seen  that  the  conditions 
T  >  Tsat(p), 

Pv  >  P> 

and 

Tv  ^  Tgat^Pv^ 


are  satisfied. 


Equation  12  gives  the  minimum  amount  of  superheat  required  for  a 
bubble  of  radius  R  to  grow.  Any  bubble  with  smaller  radius  or  less  super¬ 
heat  will  collapse.  In  Fig.  4,  we  have  compared  the  superheat  requirements 
for  water  and  for  liquid  sodium.  It  can  be  seen  that  liquid  sodium  requires 
more  superheat  for  the  same  bubble  radii. 


Fig.  4.  Growth  and  Collapse  Regions  for  Sodium  and  Water-vapor  Bubbles 


P 


Fig.  5 


E .  Growing  Process  on  Saturation. 

Curve 

We  are  considering  the  growth 
of  a  vapor  bubble  in  a  liquid  that  is 
at  uniform  temperature  and  pressure. 
It  is  therefore  assumed  that  the  tem¬ 
perature  T  and  pressure  p  remain 
constant  during  the  process. 

When  a  bubble  grows,  its 
radius  R  increases,  and  the  pres¬ 
sure  pv  and  temperature  Tv  decrease. 
In  addition,  we  know  that  the  thermo¬ 
dynamic  state  of  vapor  is  very  close 
to  the  saturation  curve.  The  state  of 
the  vapor  therefore  moves  along  the 
direction  as  shown  in  Fig.  5. 


Equilibrium  State  for 
Inertia-controlled  Growth 


The  main  task  of  bubble  dynamicist  is  to  determine  the  state  of  the 
vapor  as  a  function  of  time.  This  requires  simultaneous  solution  of  the  con¬ 
tinuity,  momentum,  and  energy  equations  for  given  initial  conditions.  The 
problem  is  very  difficult,  as  many  times  the  initial  conditions  are  not  very 
well  known. 

Let  us  look  at  the  two  extreme  cases. 

1.  Inertia-controlled  Growth 


have 


When  the  state  of  vapor  is  very  close  to  E,  as  shown  in  Fig.  5,  we 


Pv  ~  Psat(T)< 


(13) 


and 


Tv  -  T.  (14) 

This  means  that  we  have  a  larger  pressure  difference  (pv  -  p)  and  negligible 
temperature  difference  (T  -  Tv)  acting  on  a  spherical  interface.  Therefore, 
heac-transfer  rates  to  vapor  can  be  neglected.  The  growth  of  a  bubble  is 
mainly  controlled  by  inertia  force.  The  solution  of  the  momentum  equation 
alone  may  be  sufficient  to  determine  the  growth  rate. 

P  2.  Heat-transfer- 

f  controlled  Growth 


When  the  state  of  vapor 
is  very  close  to  G,  as  shown  in 
Fig.  6,  we  have 

Pv  -  P  (15) 

and 

Tv  “  T,at(p).  <16> 

In  this  case,  we  have 
a  smaller  pressure  difference 
(pv  -  p)  and  larger  temperature 
difference  (T  -  Tv).  The  heat- 
transfer  rate  is  very  high  and  can¬ 
not  be  neglected.  However,  the  inertia  forces,  being  small,  can  be  neglected. 
The  solution  of  the  energy  equation  alone  may  be  sufficient  to  determine  the 
growth  rate.  The  process  in  this  domain  is  heat-transfer  controlled. 


Fig.  6.  Equilibrium  State  for  Heat- 
transfer-controlled  Growth 
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III.  DYNAMICS  OF  A  STATIONARY  BUBBLE 

In  the  last  25  years,  considerable  work  has  been  done  to  determine  the 
growth  rate  (or  collapse  rate)  of  a  stationary  bubble.  Most  of  the  theoretical 
analyses  are  based  on  assuming  a  spherical  bubble  and  uniform  vapor  prop¬ 
erties.  Even  with  these  assumptions,  the  problem  remains  formidable.  It 
now  requires  a  solution  of  one- dimensional  conservation  equations  for  a  sur¬ 
rounding  liquid  with  suitable  initial  and  interface  boundary  conditions.  In  this 
section,  we  briefly  review  the  results  of  some  of  the  theoretical  analyses  after 
describing  the  appropriate  governing  equations. 

A.  Governing  Equations 


In  presenting  the  governing  equations,  we  have  assumed  a  spherical 
symmetry  and  considered  the  bulk- viscosity  coefficient  to  be  zero. 

1.  Conservation  Equations  for  a  Surrounding  Liquid  Phase 


a.  Continuity 


bt 


1  A 

r*2  5  r 


(r2pur)  =  0. 


For  an  incompressible  fluid,  Eq.  17  simplifies  to 
r2ur  =  constant  =  R2ur, 


(17) 


(18) 


or 


ur  =  R2uR/r2. 


(18a) 


b.  Momentum  Equation 


bur 

bt 


+  u  r 


bur 

br 


I  4.  1  (4V  X  fi\  &  f6Ur  ■  ZM 
p  br  p  V  3  /  br\  br  r  / 


+  b. 


(19) 


Here,  b  is  the  body  forces,  and  3  is  the  second  coefficient  of  viscosity.  If 
the  surrounding  fluid  is  incompressible, 


bu„  2u^ 

— ?  +  — i 

br  r 


0 


from  Eqs.  18,  and  Eq.  19  simplifies  to 


bu^ 

bt 


+  Ui 


bur 

br 


(20) 


15£  +  b. 

P  br 


(21) 
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If  body  forces  are  not  present,  Eq.  21  can  be  integrated  with 
the  aid  of  Eq.  18a  to  obtain 

2uRR  +  RuR  j  UR  =  l(pR  -  p»),  (22) 

where  •  represents  differentiation  with  respect  to  time. 

If  we  further  assume  that  the  rate  of  evaporation  is  negligible, 
then  the  conservation  of  mass  at  the  interface  leads  to  the  relation  uR  =  R 
(Eq.  30  below).  The  momentum  equation  (Eq.  22)  now  simplifies  to 

RR  +  |r2  =  -(pR  -  poo),  (23) 


which  is  Rayleigh's  equation.  With  the  aid  of  momentum  balance  at  the  inter¬ 
face  (see  Eq.  33  below),  the  motion  of  bubble  surface  can  be  expressed  in 
terms  of  vapor  pressure  pv.  Thus, 


RR  +  \  R2  = 


1 

P 


Pv 


c.  Energy  Equation 


2ur\  WSur 
r  /  3  \  br 


2  bT\ 
r  br / 


+  Pq'", 


(24) 


(25) 


where  q'n  is  the  heat  source  and  e  is  the  specific  internal  energy.  If  we 
assume  no  heat  source,  no  viscous  dissipation,  negligible  kinetic  energy,  and 
an  incompressible  fluid  with  constant  specific  heat,  Eq.  25  simplifies  to 


6  T  bT  fbZ T  2  6T\ 

Ti  +Ur  #7  =  + 

where  Of  =  k/pc  is  the  thermal  diffusivity. 

2.  Interface  Relations  and  Boundary  Conditions  at  r  =  R(t) 


(26) 


a.  Conservation  of  Mass.  If  we  consider  the  total  mass  balance 
of  a  bubble,  we  get 

Jt(^nR3pv)  =  4ttR2p(R  -  UR)>  (27) 

which  can  be  rewritten 


uR  =  eR  -  eR/3, 


(28) 
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where 


,  Pv  /-> 

e  =  1  -  — .  (2 

P 

If  it  is  assumed  that  pv  «  P,  then  e  =“  1  and  e  =“  0,  and  Eq.  28  simplifies  to 
uR  =  R.  (3 

Equation  30  implicitly  assumes  that  the  rate  of  evaporation  is  negligible. 

b.  Momentum  Balance.  For  a  concave  interface,  the  vapor 
pressure  pv  is  related  to  the  liquid  pressure  by 

Pv  =  PR  +  ^  -  Trr  +  PvuRv(uRv  "  ®-)  “  PUR(UR  '  »■)•  (3 


If  we  ignore  momentum  fluxes  at  the  interface, 


Pv  =  PR  + 


4  /^u  r  ur\ 

^\“57“  TAT 


With  the  further  assumption  of  negligible  viscous  stress,  Eq.  32  simplifies  to 


Pv  =  PR  + 


c.  Energy  Balance.  The  energy  balance  at  the  interface  leads  to 


4rrR‘ 


i?  s(*  "R3p'')[L  +  4-f 


bur  ur\l  PVR  duRv  pvR  dp, 


3  \  i)  r  r  /]  +  3  dt  3  p  v  dt 


In  Eq.  34,  it  is  assumed  that  the  kinetic  terms  are  negligible.  In  most  situa¬ 
tions,  the  latent  heat  term  usually  dominates  and  Eq.  34  simplifies  to 

=  — ;4(l"R’Pv)-  (35) 

br  4ttR2  dt\  / 


d.  Surface  Temperature.  The  boundary  condition  at  the  interface 


will  be 


Tr  -  Tv, 


as  it  is  assumed  that  the  properties  of  vapor  inside  the  bubble  are  uniform. 
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B.  Brief  Review  of  Past  Work 

It  can  be  seen  from  the  previous  section  that  the  bubble  growth  rate 
can  be  determined  from  the  solution  of  Eq,  24.  The  vapor  pressure  pv(t)  in 
Eq.  24  depends  on  vapor  temperature  Tv  (=  Tr),  which,  in  turn,  depends  on 
the  heat- transfer  rate  at  the  bubble  surface.  It  is,  therefore,  necessary  to 
solve  Eq.  26  along  with  Eq.  24. 

The  bubble  growth  rate  is  essentially  controlled  by  inertial,  thermal, 
and  diffusive  effects.  Here  we  are  considering  only  the  first  two  effects. 
Initially,  the  rate  of  growth  of  a  vapor  bubble  is  controlled  by  the  liquid  inertia, 
surface  tension,  and  the  pressure  difference  (pv  -  pm).  At  the  start  of  the 
expansion  of  the  original  bubble  nucleus,  the  forces  are  nearly  in  equilibrium 
and  the  bubble  growth  rate  is  slow.  As  the  bubble  starts  growing,  the  surface 
tension  force  reduces  and  the  growth  rate  is  accelerated.  Within  a  very  short 
time,  the  growth  becomes  appreciable  and  the  temperature  and  pressure  of 
vapor  drop.  The  growth  rate  now  decreases  and  is  essentially  controlled  by 
the  thermal  effects.  The  bubble  will  either  continue  growing  or  start  collaps¬ 
ing,  depending  on  whether  the  surrounding  fluid  is  superheated  or  subcooled. 
The  qualitative  description  of  this  process  is  shown  in  Figs.  7  and  8. 

When  the  inertia  forces  are  dominating,  the  changes  in  vapor  pressure 
can  be  neglected  and  the  solution  of  Eq.  24  is  sufficient  to  predict  the  bubble 
growth  rate.  The  results  are  valid  in  the  initial  stages  of  growth  and  for  low 


Fig.  7.  Variation  of  Growth  Rate  in  Superheated,  Saturated,  and  Subcooled  Regions 


Fig.  8.  Variation  of  Bubble  Radius  in  Superheated,  Saturated,  and  Subcooled  Regions 

Jakob  numbers.  On  the  other  hand,  when  the  thermal  effect  is  dominant,  the 
solution  of  Eq.  26  and  the  interface  relation  (Eq.  35)  are  sufficient  to  deter¬ 
mine  the  growth  rate.  Since  the  simultaneous  solution  of  Eqs.  24  and  26  has 
been  difficult,  most  of  the  analyses  are  restricted  either  to  inertia- controlled 
or  to  heat- transfer- controlled  growth. 

1.  Inertia-controlled  Growth 


Rayleigh3  was  one  of  the  first  researchers  in  bubble  dynamics. 

He  assumed  a  uniform  pressure  difference  and  integrated  Eq.  23  to  predict 
that  the  radius  of  a  bubble  increases  linearly  with  time;  i.e., 

R  66  t  (Rayleigh's  solution).  (37) 

This  solution  is  applicable  when  the  bubble- wall  motion  is  slow 
and  the  liquid  can  be  treated  as  incompressible.  Dergarabedian4  analyzed 
Rayleigh's  equation  accounting  for  the  effect  of  surface  tension,  but  assumed 
a  constant  value  for  vapor  pressure.  Hunter,5  Hickling  and  Plesset,6  and 
Hsieh7  extended  the  solution  of  the  momentum  equation  to  account  for  com¬ 
pressibility  effects. 

2.  Heat- diffusion- controlled  Growth 


a.  Uniform  Temperature  Field.  Plesset  and  Zwick8  derived  an 
integro- differential  equation  by  combining  the  momentum  equation  (Eq.  24) 
and  the  solution  of  the  energy  equation  (Eq.  26)  with  an  assumption  of  a  thin 


thermal- boundary  layer  near  the  bubble  wall.  They  obtained  an  approximate 
solution  by  dividing  the  growth  period  into  several  regimes  and  using  possible 
simplification  for  each  regime.  Their  solution 

R  =  cJa  V^t  (38) 

(c  =  2  /—  =  1.95) 

'  V  TT  7 

is  an  asymptotic  solution  when  a  bubble  becomes  large.  Equation  38  is  found 
to  be  in  agreement  with  the  experimental  measurements  even  near  the  lower 
limit. 


Forster  and  Zuber9  also  obtained  an  approximate  solution  of 
the  integro- differential  equation: 


R  +  R0  l°ge 


R 

Rn 


-  1 


Ri 

Rn 


:  Ja  y/oti 


(39) 


-  1 


(c  =  =  1.77). 

Here,  R0  is  the  initial  radius  in  equilibrium  with  the  surrounding  fluid,  and 
R}  is  the  radius  at  which  the  rates  of  growth  of  the  bubble  resulting  from  the 
superheat  term  and  from  the  evaporation  term  attain  an  equal  order  of  mag¬ 
nitude.  Since  R^  is  very  small,  Eq.  39  is  in  agreement  with  the  Plesset  and 
Zwick  solution  given  by  Eq.  38. 

Birkhoff  et  al.10  considered  a  similar  solution  of  Eq.  26. 

They  also  obtained  a  solution  of  the  form  given  by  Eq.  38,  with  the  coefficient 
c  being  an  implicit  function  of  the  Jakob  number: 

R  =  c(  Ja)  Ja  ^/ot.  (40) 

For  a  large  Jakob  number,  the  coefficient  c  asymptotically  approaches  the 
value  of  Plesset  and  Zwick.8 

Scriven11  extended  the  solution  of  Eq.  26  without  neglecting 
the  effect  of  density  ratio  Pv/p,  His  solution  is  of  the  same  form  as  that  of 
Plesset  and  Zwick,  except  that  the  growth  coefficient  c  is  an  implicit  function 
of  the  Jakob  number.  He  presented  graphs  and  tables  of  growth  coefficient 
for  a  range  of  Jakob  numbers  and  for  a  range  of  density  ratios. 

b.  Nonuniform  Temperature  Field.  Bankoff  and  Mikesell12  used 
the  Plesset  and  Zwick  perturbation  technique  to  obtain  the  growth  of  a  spheric 
bubble  in  initial  temperature  fields  that  were  exponential  and  ramp  functions 
of  the  Lagrangian  radial  coordinates. 
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Griffith13  considered  the  growth  of  a  hemispherical  bubble  at 
the  heated  surface.  He  obtained  the  finite- difference  solution  of  the  energy 
equation  and  applied  the  heat  flux  at  the  bubble  surface  to  determine  the 
growth  rate.  He  obtained  the  numerical  solution  for  three  Jakob  numbers 
(5.3,  1,  and  0.35)  and  for  three  values  (0,  0.67,  and  0.9)  of  dimensionless 
temperature : 

0  =  (Tsat  -  T)/(Tw  -  T). 

Savic14  used  a  model  similar  to  that  of  Griffith,  but  ignored 
the  variation  of  temperature  with  the  polar  angle.  He  applied  the  perturbation 
technique  to  solve  the  energy  equation  and  obtained  the  first  three  terms  of 
the  power  series 


R  =  1.95JaySt  +  B  —  +  D1—  , 

R  R 


(41) 


where  the  values  of  B  and  D*  depend  on  the  fluid  properties.  Savic  has  ob¬ 
tained  the  relation  for  computing  the  values  of  these  constants. 

Zuber15  assumed  that  the  heat-flux  contribution  to  bubble 
growth  is  less  in  a  nonuniform  temperature  field  than  in  a  uniform  tempera¬ 
ture  field.  He  modified  Eq.  35  by  the  amount  of  an  average  heat  flux  q^  from 
the  boiling  surface 


hPv 


dR 

dt 


uniform  field 


(42) 


He  integrated  Eq.  42,  using  the  relation  of  Fritz  and  Ende,16 


R' uniform  field 


.  AT 

k  — — - . 

yrr  at 


and  obtained  the  growth  rate 


(43) 


R 


qbAA^\ 

2k AT  /’ 


(44) 


Skinner  and  Bankoff17  analyzed  Eq.  26  for  a  general  spherically 
symmetric  initial  condition.  The  general  solution  was  applied  to  a  specific 
case  of  nucleate  boiling  on  a  heated  surface  (the  initial  temperature  distribu¬ 
tion  of  the  exponential  form),  and  the  growth  rates  plotted  for  a  range  of  sub¬ 
cooled  and  superheated  temperature  ratios. 


Han  and  Griffith18  considered  the  transient  conduction  equation 
for  a  layer  of  liquid  on  the  surface.  They  used  the  solution  of  this  equation 
and  the  heat  flux  at  the  surface  to  determine  the  growth  rate  of  a  bubble.  Their 
expression  for  bubble  radius, 
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R  =  R(t) ,  (45) 

is  somewhat  involved  for  a  simple  numerical  computation. 

Cole  and  Shulman19  compared  experimental  data  with  some 
existing  theories  for  Jakob  numbers  ranging  from  24  to  792.  They  observed 
that  the  growth  data  over  this  range  can  be  correlated  by 

R  =  2.5(Ja)3/4ySt.  (46) 

Van  Stralen20  modified  the  Plesset  and  Zwick  solution  by 
assuming  that  the  wall  superheat  is  a  function  of  time  and  that  only  part  of  the 
bubble  is  covered  by  superheated  layer.  His  expression  for  the  growth  rate  is 

R  =  2bJa(^)1/2  expHt/t,)1'2].  (47) 

where  b  is  a  correction  factor  of  value  less  than  one,  and  tj  is  the  time  at 
which  the  bubble  is  breaking  away  from  a  heating  surface. 

Mikic  and  Rohsenow21  expressed  the  nonuniformity  as  a 
function  of  waiting  time  and  applied  the  solution  of  one- dimensional  transient 
conduction  equation  to  obtain  an  expression  for  the  growth  rate: 


where  tw  is  the  waiting  period,  which  is  a  function  of  radius  of  cavity. 

3.  Inertia-  and  Heat-diffusion-controlled  Growth 

Mikic  et  al.22  extended  the  results  to  include  the  effect  of  inertia 
during  the  initial  period  of  growth  and  obtained  the  following  expression  for 
the  growth  rate: 


dR+ 

dt+ 


+  i  +  el 


l/ 2 


(49) 


R+  = 


t+  = 


tt  AT  /  be  \1/2  1 

12a\Tsat/  (Ja)5/Z 

(50) 

tt  be  ( AT)2 

12  o?Tsat  (Ja)3 

(51) 

where 
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and 


Tw  -  Tsat 

The  constant  b  is  equal  to  2/3  for  a  bubble  growing  in  an  infinite  mass  of 
liquid,  and  is  tt/7  for  a  spherical  bubble  growing  attached  to  a  surface. 

Equation  49  can  be  integrated  for  the  parameters  0  and  t^-.  For 
growth  in  a  uniformly  superheated  liquid  (tw  -*  ») ,  it  gives 


R+  =  |[(t+  +  1)3/2  -  (t+)3/2  -  1].  (53) 

Theofanous  et  al.23  extended  the  work  of  Ref.  22  by  accounting  for 
the  variation  of  vapor  density.  They  assumed  a  linear  pres  sure- temperature 
relation  for  the  saturation  curve  and  obtained  results  for  initial  to  final  vapor- 
density  ratios  up  to  20, 

Dalle  Donne  and  Ferranti24  numerically  integrated  both  the  energy 
and  the  momentum  differential  equations,  to  determine  the  growth  of  a  sodium- 
vapor  bubble.  They  took  into  account  the  variations  of  sodium  properties 
with  temperature  and  obtained  results  for  a  range  of  pressures,  superheats, 
and  Jakob  numbers.  Their  results  indicate  that,  even  during  the  initial  tran¬ 
sient  period,  the  heat  transfer  has  the  significant  influence  on  the  growth 
rate,  the  effects  of  all  other  factors  being  less  than  25%. 

Prosperetti25  reexamined  the  Pies  set  and  Zwick  solution  to  compare 
the  results  of  Dalle  Donne  and  Ferranti.  He  observed  that,  as  soon  as  the 
initial  effects  become  negligible  (i.e.,  when  the  radius  of  a  bubble  is  approxi¬ 
mately  10  times  the  initial  radius),  the  bubble  growth  can  be  described  by  a 
universal  function  applicable  to  any  liquid  and  any  superheat. 

C.  Growth  of  a  Sodium- vapor  Bubble 

To  determine  the  order  of  magnitude  of  a  sodium-vapor  bubble,  Eq.  53 
has  been  used  to  compute  the  growth  rate  of  a  vapor  bubble  in  a  uniform  tem¬ 
perature  field.  The  physical- property  values  of  sodium  were  obtained  from 
the  tables  and  graphs  of  Golden  and  Tokar.26  The  effects  of  pressure  and 
superheat  on  the  growth  rate  are  shown  in  Figs.  9  and  10.  Figure  11  compares 
the  growth  rates  of  water  and  sodium  vapor  bubbles.  It  can  be  seen  that  for 
the  same  values  of  superheat,  a  sodium- vapor  bubble  grows  faster  than  a 
water-vapor  bubble. 


Kig.  9.  Effect  of  Pressure  on  Growth  of  a  Sodium-vapor  Bubble. 

Conversion  factors:  1  atm  =  0.1  MPa;  °C  =  (°F  -  32)/1.8. 


Fig.  10.  Effect  of  Superheat  on  Growth  of  a  Sodium-vapor  Bubble. 

Conversion  factors:  1  atm  =  0.1  MPa;  °C  =  (°F  -  32J/1.8. 
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Fig.  11.  Comparison  of  Growth  Rates  of  Water-  and  Sodium-vapor 
Bubbles.  Conversion  factor;  °C  =  (°F  -  32) /l. 8. 


IV.  MOTION  OF  A  BUBBLE  IN  A  VISCOUS  LIQUID 

Considerable  theoretical  and  experimental  effort  has  been  focused  on 
the  dynamics  of  a  moving  bubble.  It  has  been  inferred  from  past  investiga¬ 
tions  that: 

1.  For  small  bubbles  at  low  Reynolds  numbers,  the  shape  of  a  bubble 
is  nearly  spherical. 

Z.  For  extremely  low  Reynolds  numbers  (Re  «  1),  bubbles  rise  like 
solid  spherical  particles.  Stoke' s  solution  for  a  flow  past  a  solid  sphere  is 
applicable  in  this  regime. 

3.  At  moderate  Reynolds  numbers,  the  motion  of  a  bubble  can  be 
considered  as  similar  to,  but  less  complicated  than,  the  motion  of  a  liquid 
droplet  in  a  liquid  medium.  We  can,  therefore,  apply  the  solutions  for  fluid 
spheres27-29  to  determine  the  motion  of  a  spherical  gas  bubble  in  the  low- 
Reynolds-number  regime. 
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4.  At  high  Reynolds  numbers  (Re  =*  1000),  bubbles  deform  from  a 
spherical  to  an  oblate  shape.  In  this  regime,  bubbles  vibrate  and  rise  along 
a  spiral  path. 

5.  With  a  further  increase  in  Reynolds  number,  bubbles  attain  an 
umbrella  shape  (spherical  cap)  that  is  quite  steady. 

The  dimensionless  parameters  generally  used  to  describe  the  charac¬ 
teristics  of  a  moving  bubble  are 


M  =  gn.4/pa3, 

(54) 

Re  =  2puooR/M-  (Reynolds  number), 

(55) 

W  =  2pu^R/a  (Weber  number), 

(56) 

F  -  u^/ZgR  (Froude  number), 

(57) 

and 


E  =  4gR2p/o/  (Eotvos  number),  (58) 

where  u^  is  the  steady  rising  velocity,  also  known  as  terminal  velocity,  of  ti 
bubble. 

A.  Governing  Equations 

The  problem  of  a  spherical  bubble  rising  with  constant  velocity  u^  is 
the  same  as  that  of  a  stationary  sphere  situated  in  a  downward  flowing  liquid 
with  a  constant  velocity  Uo,  at  infinity.  The  continuity  and  the  Navier-Stokes 
equation  of  motion  for  a  liquid  with  constant  density  p  and  viscosity  \±  are: 

Continuity: 


V  •  u  =  0.  (59) 

Momentum: 


1 

(u  *  V)u  =  g  -  —  Vp  +  vV  u. 

The  corresponding  boundary  conditions  are: 
u  ^  Uoo  as  | r |  oo? 


(60) 


u  ‘  n  =  0  on  r  =  R 


(61) 

(62) 
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and 


TrQ  =  0  on  r  =  R,  (63) 

where  n  is  a  unit  normal  vector. 

B.  Brief  Review  of  Some  Theoretical  Work 


Stokes30  considered  the  problem  of  the  very  slow  motion  (creeping 
motion;  Re  «  1)  of  a  sphere.  He  neglected  the  inertial  effects  and  obtained 
the  solution  of  the  equation  of  motion  (Eq.  60),  satisfying  Eqs.  61  and  62,  and 
the  no- slip  boundary  condition  uq  =  0  at  r  =  R.  The  final  results  for  terminal 
velocity  and  viscous  drag  were 


Tr  e 


R 


3  M-Uco 

2  IT 


sin  9  , 


(64) 


_  2  R2g(p'  -  P) 
9  p 


(65) 


and 


r  -  24 

C°  -  Ri’ 


(66) 


where  p'  is  the  density  of  a  sphere. 

Hadamard27  and  Rybczynski28  considered  the  motion  of  a  liquid  droplet. 
They  also  ignored  the  inertia  terms  from  the  equation  of  motion,  but  their 
solutions  satisfy  all  three  boundary  conditions  expressed  by  Eqs.  61-63.  Their 
final  results  were 


Uo> 


2 ( p1  -  p)gR2  p  +  p' 

3  2|i  +  3jjl! 


(67) 


(ug)r 


u„ 


2  (p  +  p') 


sin  0, 


(68) 


and 


24_  2p  +  3p' 
D  "  Re  3p  +  3p'  ’ 


(69) 


where  is  the  viscosity  of  a  droplet  fluid.  When  pt'  is  very  large  (|i'  »  M. ; 
rigid  sphere),  the  results  of  Hadamard  and  Rybczynski  are  in  agreement  with 
the  Stokes  solution.  These  results  are  also  applicable  to  gas  bubbles  (p1  ~  0) 
in  the  low- Reynolds -number  regime. 
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Boussinesq29  also  considered  the  motion  of  a  droplet.  He  assumed 
that  a  thin  layer  near  the  interface  has  a  higher  viscosity.  He  then  solved  the 
equations  of  motion  to  obtain 


2  P  -  p'  2  3|i  +  3p'  +  2e/R 
9  P  gR  Z\i  +  3|R  +  2e/R  ’ 


(70) 


and 


=  24  ly.  Mu'  4  2e/R  (71) 

D  Re  3 |JL  +  3p'  +  2e/R  v  7 

For  large  drop  radii,  Eq.  71  is  in  agreement  with  the  results  of  Hadamard 
and  Rybczynski. 

Levich31  analyzed  the  motion  of  a  gas  bubble  for  Re  »  1.  He  assumed 
that  there  is  a  thin  boundary  layer  near  the  interface  and  that  the  velocity 
field  in  this  layer  is  slightly  different  from  the  velocity  field  for  inviscid  flow 
past  a  sphere.  With  these  assumptions,  he  obtained  the  solution  of  the  equation 
of  motion  that  satisfied  all  three  boundary  conditions.  His  solution  for  the 
perturbed  velocity  field  is  in  the  form  of  an  integral.  The  drag  coefficient 
and  the  corresponding  terminal  velocity  are  computed  by  integrating  the  dis¬ 
sipation  energy: 

CD  =  48/Re  (72) 

and 

uco  =  PgR2/8|i.  (73) 

It  can  be  seen  that  the  drag  coefficient  is  twice  and  the  terminal 
velocity  is  half  the  corresponding  values  of  Stokes  r  flow. 

Moore32  assumed  that  at  high  Reynolds  numbers  the  flow  field  sur¬ 
rounding  a  spherical  bubble  is  irrotational;  i.e., 

ur  =  u00(  1  -  R3/r3)cos  0  (74) 


and 


u q  =  Uoo(  1  +  R3/2r3)sin  0.  (75) 

Although  this  velocity  field  does  not  satisfy  the  zero- shear- stre s s 
boundary  condition  of  Eq.  63,  Moore  used  this  solution  to  compute  the  drag 
coefficient  on  the  ground  that  the  drag  force  arises  entirely  from  the  normal 
viscous  stresses.  He  obtained 


Uoo  =  PgR2/6p, 


(76) 
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and 


CD  -  32/Re,  (77) 

which  appears  to  be  in  fair  agreement  with  experimental  data  at  moderate 
Reynolds  numbers. 

Similar  to  Levich,  Chao33  assumed  a  thin  boundary  layer  in  the 
vicinity  of  interface  and  the  velocity  field  slightly  perturbed  from  the  potential 
(irrotational)  flow  solution.  He  included  the  effects  of  internal  circulation 
and  also  corrected  the  continuity  relationship  used  by  Levich.  His  solutions 
for  the  perturbed  velocity  are  in  the  form  of  the  first  integral  of  the  compli¬ 
mentary  error  function.  He  also  interpreted  the  stress  field  at  the  interface 
to  obtain  the  drag  coefficient 


CD  = 


2p,' 
+  — 
U- 


0.314 


1  +  4p'/|jA 
Re1/2  / 


(78) 


which,  for  (i1  ~  0,  reduces  to 
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Re 


0.314N 

y/Re)' 


(79) 


Moore34  reexamined  the  motion  of  a  spherical  gas  bubble  to  include 
the  effects  of  pressure  forces  in  the  boundary  layer  for  computation  of  the 
drag  coefficient.  For  the  surface  component  of  the  velocity,  Moore  arrived 
at  the  equation  that  is  identical  to  Chao's  boundary- layer  equation.  Their 
solutions  are  slightly  different,  as  Chao  had  an  error  in  one  of  his  boundary 
conditions.  Moore  obtained  an  expression  for  the  drag  coefficient,  using 
the  momentum  argument,  in  which  he  included  the  contribution  from  the 
boundary  layer  and  wake  to  the  dissipation  calculation.  His  expression  for 
the  drag  coefficient  is 


CD 


2,211 
(Re) 1/2 


+  0(Re'5/6) 


(80) 


In  this  section,  we  have  restricted  our  review  of  some  theoretical 
work  about  the  motion  of  a  spherical  bubble.  At  very  high  Reynolds  numbers, 
the  shape  of  a  bubble  no  longer  remains  spherical.  Although  we  have  not 
reviewed  it  here,  considerable  work  has  been  done  to  evaluate  the  motion  of 
a  distorted  bubble.  Harper35  has  published  a  comprehensive  article  that  in¬ 
cludes  some  theoretical  results  for  spheroidal  and  spherical  cap  bubbles. 


Since  our  main  interest  is  to  study  the  behavior  of  a  sodium-vapor 
bubble,  we  have  presented  in  Fig.  1Z  a  graph  of  terminal  velocity  as  a  function 
of  radius  for  a  sodium- vapor  bubble,  and  in  Fig.  13  the  drag  coefficient  as  a 
function  of  Reynolds  number. 
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V.  DYNAMICS  OF  A  MOVING  BUBBLE 

We  have  seen  from  the  previous  two  sections  that  considerable  work 
has  been  done  for  the  cases  of  (l)  a  bubble  growing  but  stationary,  and  (2)  a 
moving  bubble  of  constant  volume.  Recently,  some  attempts  have  been  made 
to  analyze  the  growth  of  a  bubble  that  is  simultaneously  in  motion. 

The  thermal-boundary-layer  thickness  of  a  bubble  in  motion  (rising) 
is  smaller  than  that  for  a  stationary  bubble.  As  the  bubble  velocity  increases, 
the  boundary-layer  thickness  decreases,  and  the  heat-transfer  and  growth 
rates  increase.  As  the  radius  of  a  bubble  increases,  the  drag  force  increases, 
and  the  velocity  decreases.  This  change  in  velocity,  in  turn,  affects  the  growth 
rate.  Therefore,  we  can  see  that  for  a  moving  bubble,  the  equations  of  motion 
and  energy  are  coupled  and  require  simultaneous  solution. 

A.  Governing  Equations 

Consider  a  spherical  coordinate  system  fixed  at  the  center  of  a  bubble 
that  is  moving  (rising)  with  a  velocity  u^  in  an  incompressible  liquid  medium 
of  constant  viscosity.  Assume  that  the  flow  is  symmetric  with  respect  to  the 
vertical  axis  passing  through  the  center  of  the  bubble.  The  governing  equations 
are: 


Continuity: 


^  ^  2  1  i 

?  Sr  (r  ur  + 


1 


r  sin  9  S0 


(uq  sin  0 )  =  0. 


(81) 


Momentum: 


tt-+(u-v)u  =  g-*~ Vp  +  vV2u. 
St  p 


Energy: 


ST  ST  u9  ST 
St  +  Ur  Sr  +  r 


=  O' 


_L  JL 

r2  Sr 


(^ 


1  _S_/  .  ST: 

r2  sin  0  S0  \s*n  ®  S9) 


The  corresponding  boundary  conditions  are: 
u  -*  u^;  T  -*  T^  as  |r  |  -*  »; 
ur  =  uR  (Eq.  85);  t  rQ  =0; 

and 


(82) 


(83) 


T  =  Tsat  at  r  =  R. 


(84) 
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The  integration  of  mass,  momentum,  and  energy  equations  at  the  bubble 
interface  provides 


Un  =  eR  -  eR/ 3 , 


(85) 


D  =  ZttR2  /01T(Trr)Rsin  0  d0 , 


(86) 


and 


R  = 


m 


2P  vL  "o  '  'R 


sin  0  d0 , 


(87) 


where  t  is  the  normal  stress,  D  is  the  drag  force,  and 
r  r 


(88) 


If  we  ignore  the  change  in  the  inertia  of  a  bubble,  the  drag  force  can 
be  equated  to  the  buoyancy  force;  then, 


D  =  yTTR3g(p  -  pv) 


(89) 


or 
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2  P  UrorrR2  3  \il 


(90) 


B .  Past  Work 


Ruckenstein36  used  only  the  energy  equation  to  determine  the  heat- 
transfer  rate  between  a  translating  vapor  bubble  and  the  surrounding  liquid. 
He  assumed  a  thin  boundary  layer  near  the  interface  and  a  steady  velocity 
field  corresponding  to  the  inviscid  flow  solution.  Following  the  method  of 
Levich,  he  solved  the  steady-state  energy  equation  to  obtain 


Q  =  8ATkl 


/n\1/2/  uoo  \1/2  2 

(3)  \2Rq7  R 


(91) 


for  the  overall  heat-transfer  rate,  and 

R5/4  =  Rq/4  +  1.275JaVk  (eg/CD)l/4t 


(92) 


for  the  bubble  radius.  Here,  R0  is  the  bubble  radius  at  the  time  of  detachment. 


Tokuda  et  al.37  considered  the  transient  energy  and  diffusion  equations 
to  study  the  dynamics  of  a  moving  bubble.  They  assumed  a  potential  flow 
model  for  the  velocity  field  around  the  bubble  and  used  the  coordinate  pertur¬ 
bation  technique  to  develop  two  asymptotic  (small-time  and  large-time)  ex¬ 
pressions  for  the  growth  rate.  They  integrated  these  expressions  numerically 
to  obtain  the  radius  as  a  function  of  time.  Their  analysis  confirmed  that  a 
translative  bubble  motion  causes  a  significant  increase  in  the  growth  rates 
over  those  predicted  for  stationary  growing  bubbles. 

Chao38  made  an  analysis  for  the  transient  response  behavior  of  thermal 
or  concentration  boundary  layers  for  a  liquid  droplet  moving  at  constant  ve¬ 
locity.  He  also  assumed  a  thin  boundary  layer  at  the  interface  and  a  potential 
flow  model  for  the  velocity  field,  but  included  the  effects  of  internal  circula¬ 
tion.  He  used  similarity  transformations  and  obtained  a  solution  of  the 
transient-energy  (or  mas  s  -  concentration)  equation  in  the  form  of  a  compli¬ 
mentary  error  function.  He  observed  that  the  overall  heat-transfer  rate  is  a 
function  of  the  parameter  Ucot/R,  and  that  when  this  parameter  attains  a  value 
of  unity,  the  transient  solution  is  close  (within  1.5%)  to  the  steady- state  (large¬ 
time)  solution. 

Ruckenstein  and  Davis39  analyzed  the  transient-energy  equation  to 
determine  the  effects  of  bubble  translation  on  its  growth  rate.  They  also 
assumed  a  thin  boundary  layer  surrounding  the  bubble  surface.  They  assumed 
the  velocity  field  to  be  a  sum  of  the  potential  flow  field  and  the  flow  field  due 
to  pure  radial  growth.  They  solved  the  energy  equation,  using  similarity 
transformations.  Their  expressions  for  the  growth  rate  and  bubble  radius  are 
rather  involved  integral  equations.  They  also  presented  a  quasi-steady-state 
approximation  for  the  growth  rate  that  can  be  integrated  numerically,  provided 
the  translational  velocity  can  be  specified  as  a  function  of  radius.  They 
as  sumed 

■ (£  *  «*r  (,3) 

to  compare  their  quasi- steady -state  approximation  with  the  results  of 
Florschuetz  et  al.40 

Pinto  and  Davis41  used  Basset’s  one-dimensional  equation  of  motion 

^{pv(Pr3)u  +  ip  (f  t^R3)  u]  =  (p  -  pv)fiTR3g  -  CDnR2P7 

along  with  the  Plesset  and  Zwick  equation8  to  determine  the  velocity  of  a 
growing  bubble.  The  second  term  in  Eq.  94  represents  the  momentum  of 
liquid  accelerated  with  the  rising  bubble.  Pinto  and  Davis  neglected  the  terms 
containing  pv  in  the  equation  of  motion  and  used  the  steady-state  terminal- 
velocity/  drag-coefficient  correlations 
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>D  =  “PF 


(Levich31)  for  R  ^  0.04  cm. 


and 
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8  Pg 


D  3  1.82  a 


(Peebles  and  Garber42) 


(95) 


(96) 


to  predict  the  rising  velocity  of  a  growing  bubble.  They  obtained  an  analytical 
solution  for  smaller  bubbles  and  a  numerical  solution  for  larger  bubbles. 
Their  results  show  that,  in  the  initial  stages,  the  gravitational  force  is  the 
dominant  force  which  increases  the  rising  velocity  of  a  bubble.  But  as  the 
bubble  grows  larger,  the  drag  force  increases  and,  consequently,  the  velocity 
starts  to  decrease. 

Chang43  considered  the  growth  of  a  spherical  second  phase  (bubble  or 
droplet)  as  governed  by  simultaneous  heat-  and  mas s -transfer  limitations. 

He  showed  that  by  using  the  proper  compatibility  conditions,  one  can  apply  the 
solution  for  an  uncoupled  case  (either  heat  or  mass  transfer  limited)  to  de¬ 
termine  the  growth  rate  when  both  heat  and  mass  transfer  control  the  growth 
of  a  bubble. 


C.  Order-of-magnitude  Approach 

Aleksandrov  et  ah44  used  the  order-of-magnitude  approach  to  determine 
the  rate  of  growth  of  bubbles  that  are  rising  in  a  superheated  liquid.  Since  our 
solution  is  based  on  their  final  formula,  we  are  presenting  its  brief  derivation 
here . 


Let  <5  be  the  distance  from  the  surface  where  the  temperature  of  the 
liquid  changes  significantly,  so  that,  to  an  order  of  magnitude, 


ST  ^0°  “  ^sat  AT 

-  -n  -  =*■=  - 

dr  6  <5 


(97) 


The  changes  of  liquid  temperature  along  the  bubble  surface  occur  over  a  dis¬ 
tance  that  is  of  the  order  of  a  bubble  radius  R;  therefore, 


1  ST  AT 
r  S~0  ^  "R  ' 


(98) 


From  the  hydrodynamic  solution  of  a  rising  bubble,  the  orders  of 
magnitude  of  liquid  velocity  components  near  the  bubble  surface  are 
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and 


ur  =“  0 


.  Uco  sin  9 

n  =  -  c*  U0 

e  z 


(99) 


With  these  assumptions,  we  can  estimate  the  order  of  magnitude  of 
each  term  of  the  energy  equation 
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From  the  solutions  of  the  energy  equation  for  a  growing  but  stationary 
bubble,  the  boundary-layer  thickness  6  is  very  small  cornpared  to  the  radius 
of  a  bubble;  i.e. , 


6  «  R. 


(102) 


Retaining  only  the  terms  of  the  order  of  1/62,  we  obtain 

1/2 


ST  AT 
Sr  ~  6 


AT^  —  +-^ 
\at  ffR 


(103) 


We  can  therefore  assume  that  the  average  temperature  gradient  over  the 
bubble  surface  must  be  of  the  form 


3T 

"Sr 


R 


/  m  a  bUco 

^T~  -  Tsat\^  +  q,r 


1/2 


(104) 


where  a  and  b  are  some  numerical  coefficients.  Combining  Eqs.  87  and  104, 
we  obtain 
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(105) 


For  a  stationary  bubble  growing  in  a  superheated  liquid,  Plesset  and  Zwick8 
obtained  the  asymptotic  solution  of  the  energy  equation.  Comparing  their 
solution 


dR  _  k(Tco  ~  Tgat)/  3  V 

dt  pvL  VnQft/ 


1/2 


(106) 
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with  Eq.  10  5,  we  find  that 


Levich31  solved  the  steady-state  diffusion  equation  for  a  droplet  moving 
at  uniform  velocity.  Ruckenstein  used  the  solution  of  Levich  for  heat  transfer 
between  a  moving  bubble  and  the  surrounding  liquid.  Ruckenstein' s  final  ex¬ 
pression  for  the  average  temperature  gradient  (or  concentration  gradient  for 
the  diffusion  equation)  is 


BT 

3r 


R 


-  (T 


sat 


>( 


JL 

3tt  oR 


(108) 


On  comparing  Eq.  108  with  Eq.  104,  we  obtain 

b  =  rL  (109) 

3tt 

Equation  10  5  for  the  rate  of  growth  of  a  moving  bubble  can  now  be  expressed 
as 


dR  k(Too  ^satV  3  ^  Y^2 

dt  ”  pvL  Vnat  3vaR/ 


(no) 


Although  this  analysis  provides  a  simple  relation  describing  the  influ¬ 
ence  of  velocity  on  bubble  growth  rate,  it  does  not  consider  the  effect  of  growth 
rate  on  the  velocity  of  a  bubble. 

D.  Present  Work 


1 .  Stationary  Medium 


From  the  review  presented  in  the  previous  section,  we  can  see  that 
most  investigators  have  attempted  to  solve  the  energy  equation  for  a  bubble 
that  is  rising  at  uniform  velocity  in  a  uniformly  superheated  liquid.  Pinto  and 
Davis,41  on  the  other  hand,  considered  the  equation  of  motion  to  determine  the 
effect  of  growth  rate  on  the  velocity  of  a  bubble,  but  they  neglected  the  effect 
of  velocity  on  growth  rate.  As  far  as  we  know,  no  attempts  have  been  made  to 
consider  both  the  effect  of  velocity  on  growth  rate  and  the  effect  of  growth  rate 
on  velocity.  This,  of  course,  requires  simultaneous  solution  of  coupled 
equations . 


From  the  order-of-magnitude  analysis  of  Aleksandrov  et  al. , 
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we  note  that 
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dR  =  kAT/  3  +  Zu™_\  '  (111) 

dt  pvL\TTCk't  3ttRcv/ 


or 


dR 

dt 


2  ucA1/2 
3^  ~R/ 


f(Ja;  t;  u;  R). 


(112) 


If  we  consider  Basset's  one-dimensional  equation45  of  motion  (Eq.  94)  and 
ignore  the  terms  containing  pv,  we  obtain 


du 

dt 
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3r  U 

4CD  R 


3u  dR 
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Equations  112  and  113  are  a  set  of  coupled  ordinary  differential  equations  for 
which  numerical  solutions  can  be  obtained. 


We  have  used  a  fourth-order  Runge-Kutta  integration  scheme  to 
solve  Eqs .  112  and  113.  For  the  drag  coefficient  in  Eq.  113,  we  used 
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for 


R  <  0 .07  cm, 


and 


C  =  (Peebles  and  Garber42)  (115) 

D  3  1.82a  V  ' 

for  R  ^  0.07  cm.  The  following  relations  were  used  to  obtain  the  initial 
values : 


tmit  =  0.001  s, 


/3Q,t  initV72 


Rinit  =  2JaV~  tt 
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The  initial  time  being  very  small,  the  starting  values  of  radius 
and  velocity  did  not  affect  the  final  solution.  This  can  be  seen  from  Table  I, 
in  which  we  have  presented  the  results  for  two  different  initial  velocities.  To 
check  the  accuracy  of  the  solution,  we  made  computations  with  different  time 
intervals.  A  time  increment  of  10~5  s  gave  a  fairly  accurate  and  stable  solu¬ 
tion  (see  Table  II). 


TABLE  I.  Results  for  Water  Vapor  [Ja  =  5.0;  p  =  1  atm  (0.1  MPa)] 
with  Two  Values  of  Initial  Velocity 


Radius,  cm 

Velocity, 

cm/s 

Time, 

s 

dR 

uinit  "  dt 

uinit  =  0 

dR 

uinit  -  dt 

uinit  =  0 

0.005 

0.0292 

0.0288 

4.305 

3.787 

0.1 

0.1893 

0.1889 

25.52 

25.50 

0.9 

0.7500 

0.7498 

16.99 

16.99 

TABLE  II.  Results  for  Water  Vapor  [Ja  =  5.0;  p  -  1  atm  (0.1  MPa)] 
with  Two  Values  of  Time  Increment 

Time, 

Radius , 

cm 

Velocity,  cm/ s 

s 

At  =  1 0  ”5  s  At  = 

0.5  x  10‘5  s  At 

=  1 0 “5  s  At  = 

0.5  x  lO-5  s 

0.005 

0.0292 

0.0291 

4.305 

4.301 

0.1 

0.1893 

0.1893 

25.52 

25.52 

0.9 

0.7500 

0.74999 

16.99 

16.99 

We  have  compared  our  results  with  experimental  measurements 
of  Florschuetz  et  al.40  Figures  14-16  show  that  our  predictions  are  in  close 
agreement  with  their  measurements. 

Results  have  been  obtained  for  water  and  sodium  for  a  range  of 
Jakob  numbers.  Figures  17-21  show  the  effects  of  Jakob  number  on  velocity 
and  on  growth  rates.  The  growth  rate  and  velocities  of  sodium-vapor  bubble 
are  compared  (Figs.  21  and  22)  with  the  growth  rates  and  velocities  of  water- 
vapor  bubble.  It  can  be  seen  that  the  rising  velocity  has  a  greater  influence  on 
the  growth  rate  of  a  water-vapor  bubble  than  that  of  a  sodium-vapor  bubble. 

This  is  because,  for  the  same  values  of  Jakob  number,  the  radius  of  a  sodium- 
vapor  bubble  is  larger  than  that  of  a  water-vapor  bubble.  This  produces  a 
larger  drag  force,  lowers  the  velocity,  and,  consequently,  reduces  the  effect  of 
velocity  on  the  growth  rate.  Since  the  velocity  is  lower  for  higher  Jakob  num¬ 
bers,  the  influence  of  velocity  on  the  growth  rate  decreases  with  the  increase  in 
Jakob  numbers.  This  also  can  be  seen  from  our  results  presented  in  Figs.  17- 
20  and  22.  Figure  23  shows  the  effect  of  velocity  on  the  growth  of  a  sodium- 
vapor  bubble. 


rschuetz  et  al.  for  a  Water-vapor  Bubble  of  Florschuetz  et  al.  for  an  Ethanol-vapor  Bu] 
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Fig.  21.  Effect  of  Jakob  Number  on  Rising  Velocity  of  a  Bubble  [p  =  1  atm  (0.1  MPa)] 
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2.  Moving  Medium 


When  a  bubble  is  growing  in  a  moving  medium,  the  surrounding 
pressure  p  may  vary  significantly  along  the  path  of  a  rising  bubble.  Although 
the  surrounding  pressure  does  not  explicitly  enter  the  equation  for  growth  rate, 
it  influences  the  growth  rate  by  varying  the  vapor  density. 

The  variation  of  surrounding  pressure  differs  from  case  to  case 
because  it  depends  on  the  geometry  of  a  channel,  flow  rate,' fluid  properties, 
etc.  It  is  difficult  to  obtain  a  generalized  solution  that  is  applicable  to  all  flow 
situations;  however,  we  can  present  a  set  of  ordinary  differential  equations 
that  can  be  solved  for  a  given  flow  and  geometry  situation. 

Consider  a  bubble  that  is  rising  with  a  velocity  u  in  a  moving 
medium  of  velocity  u^.  Let  the  surrounding  pressure  pm  be  a  known  function 
of  the  axial  distance  (height)  Z: 

Poo  =  Pco(Z)-  (119^ 

If  we  assume  that  vapor  is  saturated,  we  can  evaluate  vapor  density  pv  at  p 
from  the  equation  of  state;  hence, 

Pv  =  pv(Pj-  (120) 

Equations  112  and  113  are  still  applicable  for  a  moving  medium,  except  that: 

(a)  We  now  replace  u  with  the  relative  velocity 

urel  =  u  "  V  (121) 

(b)  Our  assumption  of  uniform  Jakob  number  is  no  longer  valid. 

Thus, 


dR 

dt 


_2  urel  «\1/2 
3tt  R  / 


=  f(Ja(pv);  a;  t;  urel;  R). 


(122) 


and 


du 


rel  0  3  urel  ^urel  dR 


2g  -  |CD- 


"  *  £)■ 


dt  R  R  dt 

The  bubble  velocity  u  is  related  to  axial  distance  (height)  Z  by 
dZ 


dt 


=  u. 


(123) 


(124) 


Equations  119-124  are  a  set  of  coupled  equations  that  can  be  solved 
numerically  for  a  given  situation. 
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VI.  CONCLUDING  REMARKS 

1.  In  this  report  we  have  restricted  our  attention  to  the  growth  of  a  single 
spherical  bubble. 

2.  From  the  brief  reviews  presented  here,  we  can  see  that  considerable 
work  has  been  done  to  determine  (a)  the  growth  rate  of  a  stationary 
bubble  and  (b)  the  motion  of  a  constant-volume  bubble. 

3.  In  a  real  situation,  the  growth  and  motion  (rising)  of  a  bubble  usually 
take  place  concurrently.  Due  to  the  complexity  of  interrelations,  it  is 
difficult  to  obtain  exact  solutions  of  coupled  differential  equations. 

4.  Very  little  work  has  been  done  to  determine  the  growth  of  a  moving  bubble. 

5.  Experimental  measurements  have  confirmed  that  the  growth  of  a  moving 
bubble  is  greater  than  that  of  a  stationary  bubble. 

6.  We  have  presented  (see  Sec.  V.D)  a  set  of  coupled  differential  equations 
to  determine  the  growth  and  velocity  of  a  moving  bubble.  Our  numerical 
solution  of  these  equations  is  in  agreement  with  the  experimental  data  of 
Florschuetz  et  al.40 

7.  Results  have  been  obtained  for  a  sodium-vapor  bubble  and  a  water-vapor 
bubble  for  ranges  of  pressures  and  Jakob  numbers. 

8.  Initially,  the  buoyancy  force  dominates  the  motion  of  a  bubble,  and  the 
velocity  of  a  bubble  increases.  As  the  radius  of  a  bubble  increases,  the 
drag  force  increases  and  the  velocity  starts  to  decrease. 

9-  We  find  that,  for  the  same  Jakob  number, 

m  >  m 

\  dt /  j •  \  dt /  , 

'sodium  x  'water 

and 

(velocity)sodium  <  (velocity)water- 

10.  The  effect  of  velocity  on  the  growth  of  a  sodium-vapor  bubble  is  not  very 
large  because  the  velocity  of  a  sodium-vapor  bubble  is  comparatively 
very  low. 

11.  Although  considerable  work  has  been  done  in  bubble  dynamics,  many 
phenomena  relating  to  this  subject  are  not  understood,  e.g.,  bubble- to- 
bubble  interaction,  bubble- to- surface  interaction,  the  effect  of  nonuniform 
temperature  field  on  the  growth  of  a  bubble,  and  thermodynamic  state  at 
the  start  of  nucleation.  Further  research  in  this  area  will  generate 
answers  to  some  of  these  questions. 
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